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g Abstract 

We calculate squared helicity amplitudes for non-diagonal and mixed squark pair production at 
hadron colliders, taking into account not only loop-induced QCD diagrams, but also previously 
unconsidered electroweak channels, which turn out to be dominant. Mixing effects are included 
for both top and bottom squarks. Numerical results are presented for several SUSY benchmark 
, scenarios at both the CERN LHC and the Fermilab Tevatron, including the possibilities of light 

^) ' stops or sbottoms. The latter should be easily observed at the Tevatron in associated production 

O . of stops and sbottoms for a large range of stop masses and almost mdependently of the stop mixing 

angle. Asymmetry measurements for light stops at the polarized BNL RHIC collider are also briefly 
CLi! discussed. 
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I. INTRODUCTION 



One of the most promising extensions of the Standard Model (SM) of particle physics is 
the Minimal Supersymmetric Standard Model (MSSM) , which postulates a symmetry 
between fermionic and bosonic degrees of freedom in nature and predicts the existence 
of a fermionic (bosonic) supersymmetric (SUSY) partner for each bosonic (fermionic) SM 
particle. Since SUSY and SM particles contribute to the quadratic divergence of the mass of 
the Higgs boson with equal strength, but opposite sign, the MSSM can, inter alia, stabilize 
the electroweak energy scale with respect to the Planck scale and thus propose a solution to 
the hierarchy problem. 

Unfortunately, SUSY particles still remain to be discovered. Their masses must therefore 
be considerably larger than those of the corresponding SM particles, and the symmetry is 
bound to be broken. In order to remain a viable solution to the hierarchy problem, SUSY 
can, however, only be broken via soft mass terms in the Lagrangian, with the consequence 
that the SUSY particle masses should lie in the TeV range and thus within the discovery 
reach of current and future hadron colliders such as the Tevatron and the LHC 

Due to their strong coupling, squarks should be abundantly produced at hadron collid- 
ers. In addition, phase space favors the production of the lighter of the two squark mass 
eigenstates of identical flavor, which are superpositions of the left- and right-handed helicity 
eigenstates. Since the off-diagonal elements of the two-dimensional squark mixing matrix 
are proportional to the mass of the corresponding SM quark, squark mixing is particularly 
important for third generation squarks. As a consequence, hadropro duct ion of light stop 

and, more recently, also light sbottom pairs has received particular theoretical and 
experimental attention. Mixing effects have also been analyzed recently for slepton, and in 
particular stau, hadroproduction 

As has been observed previously. Quantum Chromodynamics (QCD) alone does not allow 
for non-diagonal (light plus heavy) squark hadroproduction at tree-level, i.e. 0{al) in the 
strong coupling constant ag- Non-diagonal squark production rather requires the presence 
of at least one squark-mixing vertex, such as the four-squark interaction, present, e.g., in 
final-state rescattering of gluon-induced diagonal squark production. The corresponding 
finite one-loop diagrams at O(a^) have been evaluated in a decoupled gluino scenario with 
top squark loops only and found to be suppressed with respect to diagonal light/heavy stop 
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pair production at the Tevatron/LHC by three to six orders of magnitude jSj. 

In this Paper, we investigate the importance of electroweak channels for non-diagonal and 
mixed squark pair production at hadron colliders. Naively, one expects these cross sections, 
which are of 0(a;^) in the fine structure constant a, to be smaller than the diagonal strong 
channels by about two orders of magnitude. For non-diagonal squark production, the in- 
terplay between loop suppression in QCD and coupling suppression in the electroweak case 
merits a detailed investigation. In the presence of the mixing of bottom squarks, their loop 
contributions must also be taken into account. Mixed top and bottom squark production is 
possible at tree-level only through an s-channel exchange of a charged VT^-boson. Observa- 
tion of this channel may allow for interesting conclusions on the supersymmetric version of 
the CKM matrix. 

Our analytical calculations are presented in Sec.|nJ where we put in the additional effort to 
calculate squared amplitudes for definite initial parton helicities to exhibit clearly the mixing 
angle dependence in the left- and right-handed components of the electroweak currents and 
to allow for future applications to polarized hadron collisions. To be complete, we present 
our results for diagonal, non-diagonal, and mixed top and bottom squark hadroproduction 
with quark and gluon initial states up to 0(a;^) and 0(q;^), respectively. 

In Sec. IIIH we choose several typical SUSY mass spectra arising from different SUSY 
breaking scenarios, including those allowing for light stops and sbottoms, and apply them 
to the Tevatron and the LHC Cross sections are presented as functions of the squark masses, 
mixing angles, and general SUSY breaking mass parameters. We also briefly discuss possible 
asymmetry measurements for light stops at the polarized RHIC collider. Our conclusions 
are given in Sec. IIVI and squark mixing is discussed in App. |X| 

It is not the aim of this work to present a full signal-to-background analysis of non- 
diagonal and mixed squark production, as this is best done within the experimental collab- 
orations and using full detector simulations. Rather, our analytical results lend themselves 
easily to irnplementation in general purpose Monte Carlo programs such as PYTHIA or 
HERWIG jlf, which are traditionally employed for these kinds of simulations. 
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II. ANALYTICAL RESULTS 



In the following, we present the leading contributions in the strong (a^) and electromag- 
netic (a) coupling constant to the color-averaged cross sections d&h^^h^ for definite helicities 
ha and hb of the initial partons a and b. We define the square of the weak coupling constant 
dw ~ ^ I si^^ in terms of the electromagnetic fine structure constant a = (? j (47r) and 
the squared sine of the electroweak mixing angle xw = siia^Ow- The coupling strengths of 
left- and right-handed (s) quarks to the neutral electroweak current are then given by 

Lg = 2T^ — 2eqXw and Rq = —2eqXw, (1) 

where the weak isospin quantum numbers are Tjj = ±1/2 for left-handed and = for 
right-handed up- and down-type (s) quarks, and their fractional electromagnetic charges are 
denoted by e^. 

In general SUSY breaking models, where the squark interaction eigenstates are not iden- 
tical to the respective mass eigenstates, the coupling strengths Lg and Rq must be multiplied 
by SjiS*i and Sj2S*2, respectively, where i,j G {1, 2} label the squark mass eigenstates (con- 
ventionally rriq-^ < rriq^) and S represents the unitary matrix diagonalizing the squark mass 
matrix (see App. In the following, the dependence on the squark mixing angle 6q will, 
however, be presented explicitly. 

Our results for the strong or electroweak 2^2 scattering processes 

ahAPa)bh,{pb) qi{pi)qf* {P2) (2) 

with a, 6 = q,q,g will be expressed in terms of the squark masses m.w , the conventional 
Mandelstam variables, 

s= iPa+ Pbf , t={Pa- Pif , and u={pa- P2Y, (3) 

and the masses of the neutral and charged electroweak gauge bosons mz and m^y. Un- 
polarized cross sections, averaged over initial spins, can then easily be derived from the 
expression 

= 1 , (4) 

while single and double polarized cross sections, including the same average factor for initial 
spins, are given by 

dAcTL = (5) 
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FIG. 1: Tree-level QCD Feynman diagrams for diagonal squark production at hadron colliders. 
The upper right diagram is absent for top squarks due to the negligible top quark density in the 
proton. 

and 

dAfTLL = ^ (6) 

so that the the single and double longitudinal spin asymmetries become 

Al = —p— and All = — — • (7) 
da da 



A. Diagonal top and bottom squark production 

For diagonal top and bottom squark production, we consider only the dominant tree- 
level strong coupling channels shown in Fig. ^ For next-to-leading order QCD corrections 
to diagonal top squark production see Ref . j3| . Other exchanges are significantly suppressed 
by smaller electroweak or Yukawa couplings. For the two upper diagrams, initiated by 



5 



quarks, we find 
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(8) 



where the upper sign holds for hi and the lower sign for 62 production. Stops are produced 
only through s-channel gluon exchange due to negligible top quark parton density functions 
(PDFs) in the proton. Even for sbottom production, t-channel gluino contributions are 
suppressed by small bottom PDFs and the heavy gluino propagator, t^^ = {t — m?)~\ 
where m?. is the gluino mass. In the case of no squark mixing, our results agree with the 

n 

double-polarized and unpolarized cross sections in Ref. 

For the four lower diagrams in Fig. ^ initiated by gluons, we find 

"2 r / I 



dt 



where ta 
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(9) 



u 



independently of the squark mixing angle and thus in 



direct agreement with the double-polarized and unpolarized cross sections of Ref. 



B. Non-diagonal top and bottom squark production 

Due to the non-chiral strong interaction of the gluon, non-diagonal squark pairs cannot be 
produced in QCD at tree-level with either qq or gg initial states. There are, however, quark- 
induced chiral tree-level electroweak and supersymmetric interactions, proceeding either 
through an s-channel Z-boson exchange (see Fig. |21) or a t-channel gluino or neutralino 
exchange (see Fig. IHI). As in the case of diagonal squark production, t-channel exchanges are 
present only for bottom squarks, due to negligible top quark PDFs, and even there they are 
suppressed by small bottom PDFs and the heavy gluino or neutralino propagators, where 
the latter are defined by tzl = (t — m?o)^^. m^o is the mass of the fc*'^ neutralino mass 
eigenstate, which couples with strength 

= ± \eq^xw{\-xw)Ki + - e.XH/) K\ (10) 
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FIG. 2: Dominant tree-level electroweak Feynman diagram for non-diagonal (i 7^ j) squark pro- 
duction at hadron colliders. 
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FIG. 3: Feynman diagrams for non-diagonal {i 7^ j) squark production at hadron colliders pro- 
ceeding through t-channel gluino (left) and neutralino (right) exchange. 

to left- and right-handed (s) quarks, if we neglect the relatively small bottom Yukawa cou- 
plings. The unitary matrix A^' diagonalizes the neutralino mass matrix in the photino-zino 
basis Gluino exchanges are, in contrast, enhanced by the strong coupling constant a^. 

Our analytical result up to Oi^a^) and 0{a1) for s- and t-channel exchanges, including 
interferences, is 
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FIG. 4: Subdominant loop-level QCD Feynman diagrams for non-diagonal [i ^ j) squark produc- 
tion at hadron colliders. 
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where we have summed over qiq^ + (?2Q^i final states. Note that there is no neutralino-gluino 
interference term due to color (no n-) conservation. 

Within QCD and in the limit of a decoupled gluino, the only possibility to produce 
non-diagonal squark pairs is by rescattering of diagonal squark pairs through four-squark 
vertices in the final state. The corresponding Feynman diagrams, shown in Fig. |3J have 
(one-) loop topology and are therefore suppressed by additional squark propagators and/or 
annihilations within the squark loop. The squared helicity amplitude for the production of 
non-diagonal stop pairs in gluon-gluon collisions is given by 
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FIG. 5: Dominant electroweak Feynman diagram for mixed {e.g. top and bottom) squark produc- 
tion at hadron colliders. 



and Pq. = Y 1 — 4m|. is the velocity of the i squark mass eigenstate. In the limit of 
degenerate light squarks, only top and bottom squark loops survive loop annihilations, and 
the squared helicity amplitude simplifies to 



dt 
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(14) 



23047rs4 

These expressions have been summed over tit^ + ^2^1 final states and generalize the cor- 
responding result in Ref. j^, where only top squark loops were taken into account. For 
non-diagonal sbottom production, top and bottom squark indices have to be exchanged in 
the equations above. 



C. Mixed top and bottom squark production 



As mentioned above, mixed top and bottom squark production proceeds at tree-level only 
through an s-channel exchange of a charged IV^-boson, shown in Fig. El The exchange of 
a t-channel chargino is excluded by negligible top quark PDFs in the proton, and t-channel 
gluino contributions are loop-suppressed, since they require electroweak rescattering of the 
intermediate squarks. The squared helicity amplitude is therefore 



dt 



tu 



2 2 

ti bi 



cos^ 9t cos^ 9i{l - ha){l + hb) 



(15) 



where we show explicitly the dependence on the SM and SUSY CKM matrix elements. Here, 
we have not summed over iibl + biil final states, since each of them is produced by different 



partonic luminosities of the correct weak isospin partners q and q' in initial state. This 
purely left-handed charged current cross section is easily derived from the squared s-channel 
contribution in Eq. (fTTj) by adjusting the squark masses, including the squared absolute 
values of the appropriate CKM matrix elements Vqq' and Vi^i^ and by setting 

niz —>■ mw, -Rg = 0, Lq = V2 cos^ 6w , and sin^(26'g) 4cos^ 6'f7Cos^ ^j;. (16) 

For the mixed production of the heavier squark mass eigenstates, the corresponding index 
1 has to be replaced by 2 and the squared cosine of the mixing angle by the squared sine. 
At lowest order, gluon initial states do not allow for the production of a (charged) mixed 
squark final state. 



III. NUMERICAL RESULTS 



We now turn to our numerical results for the production of diagonal, non-diagonal, and 
mixed squark pairs at hadron colliders and use the squared helicity amplitudes given in 
Sec. ini together with Eq. (jH) to calculate unpolarized partonic cross sections. The QCD 
factorization theorem then allows to compute unpolarized hadronic cross sections 

fl (•1/2 In T /-imax „ „ H(J 

= dr/ dy dt fa/A{xa,Ml) h/B{x,,Ml) — (17) 

Jm?/S J-l/21nr Jt^^^-, OX 

by convolving the relevant partonic cross section do" with universal parton densities fa/ a and 
fh/B of partons a, h in the hadrons A, B. The PDFs depend on the longitudinal momentum 
fractions of the two partons Xa,b = y/re^^ and on the unphysical factorization scales Ma^b- 
We use the most recent leading order (LO) PDF set by the CTEQ collaboration, CTEQ6L1 
jlol |. at the factorization scale Ma = Mf, = m = {rriq^ + m_{/))/2 and identify the latter with 
the renormalization scale /z in the strong coupling constant as(/i). The QCD scale parameter 
A in the CTEQ6L1 fit for Uf = 5 active quark fiavors is 165 MeV. 

For the masses and widths of the electroweak gauge bosons, we use the current values 
of mz = 91.1876 GeV, mw = 80.425 GeV, Tz = 2.4952 GeV, and Tw = 2.124 GeV. The 
squared sine of the electroweak mixing angle 

sin^6'iy = ^ — 'mw/'m% (18) 
and the electromagnetic fine structure constant 

a = m^sm'^ 6w/t^ (19) 
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can then be calculated in the improved Born approximation using the world average value 

n 

of Gf = 1.16637- 10"^ GeV-2 Fermi's coupling constant [11]. For the SM CKM matrix 



ll|, while we set the (so far unknown) SUSY 



elements, we take the central values in Ref. 
CKM matrix elements V^.^. to one. 

The physical masses of the produced squark mass eigenstates and mixing angles are calcu- 
lated using the recently updated computer program SUSPECT [3]. Its Version 2.3 includes 
now a consistent calculation of the Higgs mass, with all one-loop and the dominant two-loop 
radiative corrections, in the renormalization group equations, that link the restricted set of 
SUSY breaking parameters at the gauge coupling unification scale to the complete set of 
observable SUSY masses and mixing angles at the electroweak scale. We choose two recently 
proposed minimal super gravity (mSUGRA) points, SPS la and SPS 5, as benchmarks for 
our numerical study 3]. 

SPS la is a typical mSUGRA point with an intermediate value of tan (3 = 10 and ^ > 0. 
It has a model line attached to it, which is specified by ttiq = — = 0.4 mi/2- We vary 
the common fermion mass mi/2 from 100 GeV, where rrif^ = 177 GeV and m^,^ = 229 GeV 
lie already considerably above the current exclusion limits of 95.7 and 89 GeV \v^, to 300 
GeV for the Tevatron (Fig. and and 500 GeV for the LHC (Fig. [7j), respectively. At the 
benchmark point, mi/2 = 250 GeV, leading to relatively heavy ti and bi masses of 399 and 
521 GeV. 

With the limited Tevatron center-of-mass energy of 75 = 1.96 TeV, non-diagonal and 
mixed squark production will be difficult to discover. As one can see in Fig. IHl only diagonal 
production of the lighter top squark mass eigenstate will be visible in the full region of the 
mSUGRA parameter space shown here with the expected final integrated luminosity of 8.9 
fb~^. For diagonal sbottom production, the accessible parameter space is already reduced 
to mi/2 < 225 GeV. Non-diagonal and mixed squark production could only be discovered, 
if the common fermion mass mi/2 is not much larger than 100 GeV. As expected for a 
pp collider, the cross sections are very much dominated by qq annihilation, even for the 
diagonal channels, and gg initial states contribute at most 15% in the case of diagonal light 
stop production. 

The LHC with its much larger center-of-mass energy of = 14 TeV and design lumi- 
nosity of 300 fb~^ will, in contrast, have no problem in producing all combinations of squarks 
in sufficient numbers. The hierarchy between the strong diagonal production channels of 
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p p ^ q. q'.^, Tevatron, \Is= 1.96 TeV 
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FIG. 6: Production cross sections for top (full), bottom (dashed), and mixed top and bottom 
squarks (dotted) at the Tevatron as a function of the common fermion mass /2 in the mSUGRA 
model line SPS 
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0{a1) and the electroweak non-diagonal and mixed channels of 0{a^) is, however, clearly 
visible in Fig. [7[ the latter being about two orders of magnitude smaller. The LHC being 
a pp collider and the average x-value in the PDFs being considerably smaller, the diagonal 
channels are enhanced by the high gg luminosity, which dominates their cross sections by 
up to 93 %. Among the electroweak 0{a^) processes, mixed production of top and bottom 
squarks is favored over non-diagonal top or bottom squark production by the possibility of 
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FIG. 7: Production cross sections for top (full), bottom (dashed), and mixed top and bottom 
squarks (dotted) at the LHC as a function of the common fermion mass mi/2 in the mSUGRA 
model line SPS la Q. 



two light masses and a positive charge in the final state, which is more easily produced by 
the charged pp initial state. 

SPS 5 is a slightly different mSUGRA scenario with lower tan/? = 5, larger mi/2 = 350 
GeV, and large negative Aq = —1000 GeV, leading to heavier sbottoms of 566 and 655 
GeV, a heavy t2 of 651 GeV, but also a light ti of 259 GeV. Since non-diagonal and mixed 
squark production will be unaccessible in this case at the Tevatron, we show in Figs. |H1 and IHl 
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pp^q- q\ * LHC, "{s = 14 TeV 




[GeV] 

FIG. 8: Production cross sections for top (full), bottom (dashed), and mixed top and bottom 
squarks (dotted) at the LHC as a function of the scalar mass mo in the mSUGRA model SPS 5 
|13| with a light top squark. 

numerical results for the LHC only, varying either mo (Fig. |H)) or Aq (Fig. ^ independently 
to test the sensitivity of the cross section on the squark masses and mixing. 

Concentrating first on the mo- dependence, we see in Fig. |Hla clear hierarchy between the 
dominant pair production of the lighter stop, strong pair production of the heavier stop and 
sbottoms, charged and neutral electroweak production of final states involving at least one 
light squark, and finally charged and neutral electroweak production of the heavier squarks. 
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which may only be visible up to mo < 600 GeV. The more pronounced hiearchy in Fig. |Hlcan 
be explained by the considerable squark mass differences in SPS 5, which lead to additional 
phase space suppression for the heavier squarks. 

As mentioned in the Introduction, the case of non-diagonal squark production merits a 
more detailed investigation of the relative importance of the coupling-suppressed electroweak 
diagram in Fig.|21and the loop-suppressed QCD diagrams in Fig. HJ As mentioned in Sec. lTTl 
the t-channel gluino and neutralino exchange diagrams in Fig.Elare only present for sbottom 
production. For the mSUGRA scenario considered here, their contributions are found to be 
six to eight orders of magnitude smaller than the s-channel contribution in Fig. El While one 
naively expects the and 0{af) diagrams in Figs. |21 and 0] to contribute with roughly 

equal strength. Fig. El shows for the case of non-diagonal stop production in the SPS 5 
scenario and at the LHC that the QCD loop-contributions are smaller than the electroweak 
tree-level contribution by about one order of magnitude. This is easily explained by the 
presence of additional heavy squark propagators in the loop diagrams. Here, we consider 
not only loops involving top, but also bottom squarks, which do not cancel in Eq. (fT^ . if the 
masses of the two sbottom mass eigenstates are unequal. However, the non-diagonal elements 
in the squark mass matrices are proportional to the relevant SM quark mass and rrih <^ m^, 
so that mixing effects are less important for sbottoms than for stops. Consequently, sbottom 
loops contribute about one order of magnitude less than stop loops, as can also be seen in 
Fig. El SUSY-QCD loop diagrams involving gluino exchanges have not been calculated here, 
as they are of O^a^) and require in addition the presence of heavy top quark and gluino 
propagators in the loop. In the SPS 5 scenario, we have indeed a heavy gluino of mass rrig = 
725 GeV. The tree-level cross section in Eq. (Illj) and the sbottom loop contribution in Eq. 
(fT^ depend on the cosine of the stop mixing angle through sin^(26'j) = 4 cos^ 9f (1 — cos^ 9f), 
which is clearly visible in Fig. El In contrast, the stop loop contribution in Eq. (fT^ has a 
steeper dependence through sin^(46'(-) = 16cos^^^f (1 — cos"^ Oi) (1 — 2cos^6'i)^, which is also 
visible in Fig. El 

The apparent excess of the experimentally observed bottom cross section at the Tevatron 
over the QCD prediction has led to speculations of light sbottom contributions to this cross 
section ^, although the use of up-to-date information on the B fragmentation function may 
reduce the discrepancy to an acceptable level Taking the light sbottom scenario at 

face-value with a sbottom mixing angle of sin 6*7 = 0.38 that reduces its coupling to the Z-, 
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FIG. 9: Contributions of tree- level electroweak and loop-level QCD processes to non-diagonal stop 
production at the LHC in the mSUGRA model SPS 5 3| with a light top squark as a function of 
the cosine of the top squark mixing angle, together with their sums. 

but not the ly-boson, we show in Fig. that mixed production of light top and bottom 
squark mass eigenstates at the Tevatron is yet another promising channel to confirm or 
exclude this scenario, as the cross section is well visible over a large range of hi masses, 
almost independently of the stop mixing angle and for light {mi^ = 100 GeV) as well as for 
heavier (m^^ = 400 GeV) top squarks. 

As a possible application of the polarization dependence of our analytical results, we show 
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FIG. 10: Associated production of top and bottom squarks at the Tevatron as a function of a 
light sbottom mass jj] for m^^ = 100 GeV (dashed) and 400 GeV (full) and no (top) or maximal 
(bottom) top squark mixing. 

finally in Fig. ^2 the double-spin asymmetry of diagonal light stop production at RHIC. As 
this polarized pp collider has only a rather small center-of-mass energy of = 500 GeV, 
the observable stop mass range is obviously very limited. Already for m^-^ > 106 GeV, the 
unpolarized cross section drops below 1 fb, while stop masses below 96 GeV are most likely 



already excluded 



This leaves only a very small mass window of 10 GeV for possible 



observations. In Fig. ^2 one can clearly see the rise of the asymmetry for qq and gg initial 
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FIG. 11: Contributions of the qq (dotted) and gg (dot-dashed) initial states to the longitudinal 
double-spin asymmetry An at RHIC, together with their sum (full), for stop pair production as 
a function of the light top squark mass. The total asymmetry using GRSV-VAL (dashed) instead 

n 

of GRSV-STD (full) parton densities J^] is also shown. 

states, as the stop mass and the correlated x-value in the PDFs grows. However, as the 
two asymmetries are approximately of equal size, but opposite sign, the total observable 
asymmetry rests below the 5% level in the entire mass range shown. This is true for both 
choices of polarized parton densities, GRSV 2000 standard (STD) as well as valence (VAL) 

cross sections have been calculated in this case using 
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the GRV 98 parton density set [1 



IV. CONCLUSION 



In Summary, we have presented in this Paper the most extensive analysis to date of 
diagonal, non-diagonal, and mixed squark production at hadron colliders. Great care has 
been taken to include in all cases the dominant contributions mediated not only by strong, 
but also by electroweak interactions up to 0{a'^) and (9(q;^), respectively. Squared helicity 
amplitudes have been presented for all considered channels in analytic form, as they expose 
the left- and right-handed contributions in the electroweak channels, allow for future appli- 
cations to polarized hadron collisions, and may easily be implemented in general purpose 
Monte Carlo programs. 

Numerically, we have focused on top and bottom squark production, including mixing 
in both cases, at the Tevatron and the LHC We have emphasized the fact that associated 
light sbottom and stop production may allow for confirmation or exclusion of light sbottom 
scenarios at the Tevatron. In more traditional scenarios such as the SPS la or SPS 5 
models, non-diagonal and mixed squark production can probably only be studied at the 
LHC, where these channels may allow for additional constraints on SUSY masses, mixing 
angles, or the SUSY CKM matrix. 



Note added in proof: After presentation of our work at the GDR SUSY conference in Greno- 
ble ^1 and at the Cortona 2005 Theoretical Physics Meeting [is!], and while this Paper 
was being completed, a second publication related to the mixed top and bottom squark 
production aspects of our paper appeared Unfortunately, the authors do not present 
any cross section formulae, so that analytical comparisons are impossible. However, we have 
checked that our calculations agree with their numerical s-channel results in Tab. I. Note 
that real emission contributions such as those labeled t-channel in Tab. I, coming from gluon 
initial states, cannot be reliably calculated separately, as they partially contain initial-state 
singularities that must be absorbed in the proton PDFs, so that these contributions must 
be included in a complete next-to-leading order calculation. 
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APPENDIX A: SQUARK MIXING 

The (generally complex) soft SUSY-breaking terms Ag of the trilinear Higgs-squark- 
squark interaction and the (also generally complex) off-diagonal Higgs mass parameter n 
in the MSSM Lagrangian induce mixing of the left- and right-handed squark eigenstates 



<lL,R of the electroweak interaction into mass eigenstates qi^2- The squark mass matrix 



+ ^1 ">'q">'LR 



= I " - \ (Al) 

mgrriLR + ml 



with 



= {Tq - Cq sin^ 9w)m\ cos 2/3 + m|, (A2) 

for up — type squarks, 
'^RR = CgSin 9wmzCos2(3 + <^ (A3) 

m|, for down — type squarks, 

, cot (3 for up — type squarks 
rriLR = Ag-fi* { (A4) 
tan (3 for down — type squarks 

is diagonalized by a unitary matrix S, SJ^'^S'^ = diag(m^,m2), and has the squared mass 
eigenvalues 



,2 2 



1,2 q ■ 2 



+ \ {iTT'LL + mlji T \J {rnlL - mlj^y + AmUmLR]"^) . (A5) 



For real values of mm, the squark mixing angle Oq, < 6q < n /2, in 



S 



can be obtained from 




w,th |=S| I (A6) 

\^ - sin 6q cos 6q J 



2m„mTT, 

tan 29q = „ " , . A7 

^LL - ^RR 
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If is complex, one may first clioose a suitable phase rotation = c^'^qr to make the 
mass matrix real and then diagonalize it for ql and g^. tan/5 = Vu/vd is the (real) ratio of the 
vacuum expectation values of the two Higgs fields, which couple to the up- and down-type 
(s)quarks. The weak isospin quantum numbers for left-handed up- and down-type (s)quarks 
with hypercharge Yq = 1/3 are Tjj = {+1/2, -1/2}, whereas Yg = {4/3, -2/3} and = 
for right-handed (s)quarks, and their fractional electromagnetic charges are = + Yq/2. 
The soft SUSY-breaking mass terms for left- and right-handed squarks are ttiq and m^, rrif), 
respectively, and mz is the mass of the neutral electroweak gauge boson Z^. 
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